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Abstract. Recent progress on molecular theories of hydration of nonpolar 
solutes in aqueous solution has led to new ways of thinking about the old 
issue of free, or "available", volume in liquids. This article surveys the 
principal new results with particular attention to general issues of packing 
in liquids. 



1. Introduction 

Aqueous solutions of colloidal solutes are preeminent examples of complex 
liquids. In such settings, attention is often directed towards the issues of 
macromolecular structure, aggregation, and dynamics. However, one as- 
pect of these problems is stubbornly associated with the small size scale 
of molecules, in particular the small size of a water molecule. That prob- 
lem is the molecular understanding and description of hydration effects - 
hydrophobic effects - that stabilize membranes, micelles, folded proteins, 
and aggregates of such structures. This presentation will focus on that ba- 
sic, molecular scale issue underlying aqueous solutions of interest both to 
biophysics and material science of colloids. 

Because hydrophobic effects are so broadly discussed, it might be sur- 
prising to notice that workers on molecular theories of hydrophobic effects 
have not achieved good agreement on that molecular theory, at least when 
we proceed beyond the primitive stage of the principles of statistical me- 
chanics and specification of the intermolecular interactions involved. The 
empirical fact is that different researchers hold different opinions on the 
correctness of several available theories, each of which agrees with a set of 
experimental data at least roughly. A sampling of different perspectives is 
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available from the references [1-33]; this collection is not intended to be 
complete, however. 

Furthermore, the simplest extensions of our experimental information 
can spark new debates about our understanding of hydrophobic effects. 
A current example, the effect of pressure on hydrophobic stabilization of 
folded proteins, is discussed below. 

In consequence of such observations, the theory surveyed here was de- 
signed for maximal simplicity on the fundamentals of statistical mechanics 
and on the physical assumptions applied to the particular problems. The 
theory that was developed is relevant to the traditional packing problems 
of theories of liquids and we emphasize that connection in the discussion 
here. 

A central problem for the theory of hydrophobic effects is the old, but 
imperfectly solved, theoretical problem of finding space for a solute in liquid 
solvents. As is well recognized, liquids arc dense, disordered materials. It 
is this combination of attributes that makes these problems difficult. The 
description of available volume in hard core model liquids is central to the 
modern understanding of the van dcr Waals equation of state and is basic to 
such ordering phase transitions as the hard sphere freezing and the liquid 
crystal phase transitions [34]. In the specific motivating case considered 
here, the particular molecular structuring characteristic of liquid water is 
expected to be important. So we must preserve the fidelity of the description 
of the structure of liquid water, in addition to tackling the problem of 
packing in dense, disordered materials. 

This lecture develops a new conceptualization of this old problem and 
new techniques for predicting the fractional free volume accessible to hard 
core molecules in condensed phases. This new approach is based upon an in- 
formation theory perspective that has general applicability and was initially 
explicitly heuristic. Information was sought on a condensed medium of in- 
terest and on that basis a prediction of the fractional free volume accessible 
to a hard model solute was made. For some important cases considered so 
far, the detail of information required for accurate, interesting predictions 
has been surprisingly modest. 

2. Available Volume Statistics 

Consider the solubility of inert gases in aqueous solutions. The medium is 
liquid water and the solute is idealized as a hard object, perfectly repelling 
the center (oxygen atom) of each water molecule. For such models the 
interaction part of the chemical potential of the solute is obtained as 
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with po the probabihty that the hard solute could be inserted into the 
system without overlap of van der Waals volume of the solvent; l//?=kBT. 
This is a specialization of Widom's formula [35, 36] 



AC/ is the change in the solute-solvent interaction potential energy upon 
placement of the solute in an arbitrary position in the solvent and the aver- 
age indicated by (• • •)o is over the thermal motion of the solvent unaffected 
by the solute. The solute is a test particle for this calculation. For the hard 
core model being considered, AC/ is either zero or infinity, so the average 
sought involves a random variable with value either one or zero; the av- 
eraging collects the fraction of solute placements that would be allowed. 
If presented with a thermal configuration of a large volume of solvent, we 
might estimate these quantities by performing many trial placements of the 
solute throughout the solvent and determining the fraction of those trial 
placements that would be allowed. This estimates Vf.,.^,f,/V , the fractional 
free volume accessible to the solute. Thus, Eq. (1) is a free volume formula 
[37], exact for the model being considered. 

The operation of these formulae can be viewed alternatively: Imagine 
identifying a molecular scale volume at an arbitrary position in the liquid 
system by (1) hypothetical placement of the solute and (2) determination 
of those positions of water oxygen atoms that would be excluded due to 
solute-solvent interactions. We will call this volume the observation volume. 
With such a molecular scale volume defined we could keep track, say during 
a simulation calculation, of the probabilities p„ that n = 0, 1, . . . oxygen 
atom occupants are observed. As the notation suggests, po is the probability 
that no occupants are observed in the molecular volume. 

Our strategy for predicting of po will be to model the distribution p„ 
and to extract the extreme value po. This is a primitive approach to theo- 
ries of and solubilities of inert gases water. Both more and less subtle 
theoretical works on these topics have been long available. "Less subtle" 
here means simulation calculations, techniques more straightforwardly use- 
ful than many "more subtle" approaches. The "more subtle" means here 
that further statistical quantities have been introduced for spherical solutes 
with the intention that they might facilitate more expansive approximate 
theories. These include 



exp{-/3A//} = (exp{-/?AC/})o . 



(2) 



di{X) = - 



dpojX) 
dX 
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A is the radius of a center-to-center exclusion sphere. These quantities have 
been useful in suggesting physical theories because they have interpreta- 
tions that arc appreciated physically. di{X), is the distribution function of 
distances A from an arbitrary point in the liquid to the nearest solvent cen- 
ter. 4TrX^pG{X) [Eq. (4)] is, in view of Eq. (1), the derivative with respect to 
exclusion radius of the hydration free energy due to intermolecular interac- 
tions, in thermal energy units k^T. It gives the compressive force exerted 
by the solvent on the hard spherical solute. In addition, AirX^ pG{X)d\ is 
the expected number of solvent centers in a shell of radius A and width dX 
outside a hard sphere that excludes solvent centers from a ball of radius A. 

3. Simulation Results for Liquid Water 

Some of the simulation work has determined the quantities of Eqs. (3) and 
(4) for molecular liquids represented realistically at the current state-of- 
the-art [14, 16]. Thus we know that di{X) for liquid water and for liquid n- 
hexane are both unimodal with maximum displaced by a distance less than 
O.lA; the maximum occurs at slightly smaller distances for liquid n-hcxane 
than for liquid water [14, 16]. This difference in the most probable cavity 
size between liquid n-hexane and liquid water is not large. The difference in 
the most probable cavity size between liquid water and a reference random 
medium, with sites of the same radius and distributed randomly at the 
water density, is greater than the difference between n-hexane and water. 
So the differences observed between the two molecular liquids considered 
are not purely reflections of molecular size and density. This comparison 
addresses the idea that the low solubility of inert gases in liquid water 
might be due to the small size of the water molecule and the possibility 
that "interstitial" cavities would, on this basis, be smaller in water than 
in the organic liquid [3, 38]. The fact that these differences between water 
and n-hexane are slight is associated with the fact that the basic units 
considered in n-hexane are the methyl and methylene groups. These are 
not so different in size from a water molecule. It should be noted also that 
on a packing fraction basis, typical organic liquids are denser than liquid 
water [6]. 

The notable distinction between the results for liquid water and liquid 

n-hcxanc is that the distribution di{X) is narrower for liquid water. This 
suggests that the liquid water phase is less flexible than the liquid n-hexane 
phase in opening cavities of substantial size. 

Simulation calculations have also produced G(A) for 0< A <3.0A, ap- 
proximately. This size range covers the simplest atomic solutes He and Ne 
but not much more. This does, however, permit comparison between water 
and organic solvents, and it permits comparison of available theories with 
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the simulation data. Thus for the range 2.0A< A <3.0A, G(A) for liquid 
water is approximately two-times larger than for n-hexane. Water exerts a 
higher compressive force on the surface of an inert solute than do typical 
organic liquids; water squeezes-out hydrophobic solutes [39]. 

The checking of theories against the available simulation data G{X) for 
water has also been revealing [14, 16]. We now know that the predictions 
of the scaled particle model [40, 41] are significantly below the numerically 
exact results for G(A). The Pratt-Chandler (PC) integral equation theory 
[42-44] predicts results for G(A) that are significantly too large. The more 
pragmatic revised scaled particle model due to Stillinger [2] typically pre- 
dicts G{X) between those two theories and with some empiricism about 
interpolation junctions can describe the available simulation data satisfac- 
torily [14, 16]. For sizes A 3>3.0A the available computer simulation data 
are less extensive, the theories less convincing, and the checking has been 
pursued less vigorously. See, however, the recent results of Reference [45]. 




Figure 1, Maycr-MontroU expansion for the insertion probability po. The standard 
notation foUows [50, 51], e.g., the sohd bonds indicate Mayer f-functions. The second hne 
gives the evaluations for the diagrams shown in the case of a hard core solute, n is the 
number of solvent centers in the observation volume. 



4. Information Model 

We now return to the theoretical program of predicting po. What are the 

standard theoretical tools for this? The most immediate guiding theory is 
the 'inclusion-exclusion' development [46] of Eq. (2) [47, 48, 49]: 

PO = 1+ E M- dvi dV2--- dVmP^^\v^,V2,...,rm) , (5) 
m=l ' 

where the m-body joint density for solvent centers. This is depicted 

in a standard way in Figure 1. These are standard combinatorial results, 
frequently seen in forms such as [46] 

PO = l-(n)o+ E M-M^-l)---("-"^ + l))o • (6) 

2 7TI. 
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- <n>o + (l/2){<n(n-l)>Q- <n>o^} + 



Figure 2. 'Virial' expansion, notation as in Figure 1. Virial typically connotes a 
low density expansion but here we are assuming full knowledge of the medium corre- 
lation functions in the absence of the solute; <5p^^'(l,2) = p^^^(l,2) - p(^'(l)p(^'(2). 
Thus from the perspective of a density expansion this organization shuffles the con- 
tributions to the virial coefficients. In fact, here contributions are ordered according 
to the number of bonds attached to the root point. Successive contributions have a 
structure that may be derived from that of familiar cumulants with the formal replace- 
ment {n'°)o {n{n — 1) ■ • • (n — fc + l))o; then po can be formally expressed as {e~")o- 
See [52, 53, 54]. Table 1 gives formulae for contributions through 5th order. 



Here the random variable n is the number of solvent centers within the 
observation volume and, e.g., (n)o is the expected number of centers within 
the observation volume. Several important points can be made from Fig- 
ure 1. The first is that the primordial available volume model is obtained 
from the first two terms shown (3Aii ~ — ln[l — (?t-)o]. The second point is 
direct and basic: po is naturally expressed in terms of occupancy moments, 
indeed binomial moments here. The sum truncates sharply for cases where 
a finite maximum number of particles can be present in the observation 
volume. The sum can be of practical value in the case where only the first 
nontrivial term is retained. For large solute volumes or solvent densities, 
this sum is not directly useful; we seek a way to exploit the same informa- 
tion but in a more broadly useful form. 

The next most immediate theoretical guidance comes from the 'virial' 
expansion depicted in Figure 2. This can be considered a resummation 
of the series Figure 1 and is better in the sense that such a truncation 
cannot produce a negative probability as truncation of Figure 1 can do. 
However, when that trouble is avoided Figure 2 can still be less compact and 
truncations can be less accurate. For example, single term approximation 
to Figure 1 can be satisfactory and then Figure 2 with a single term is likely 
to be less so. 

Ultimately these considerations avoid the issue that we have only limited 
information and we want to make the best prediction of po that we can. 
When the problem is stated this way what to do next is clear: we model 
the probabilities Pn on an information theory basis. We consider a relative 
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TABLE 1. Successive contributions to the series /3A/Lt = aj of 
Figure 2. (") = n{n — \) • • • in — j -\- is the binomial coefficient. 



j 

1 +((i))o 

2 -((2))o + ((^)>o72 

3 +((3))o-((i)>o((^)>o + ((^)>oV3 

4 -(G))o + ((i)>o((3))o + ((2))o72 - ((T))o7(2)>o + ((")>oV4 

5 +((5))o-((T)>o((:))» + (("))..7(3)>.. 

-(G))o((3))<. + (C;))o(G))o^ - ((';))o7G))» + (('i))oV5 



or cross information entropy [55], 

riiiPn}) = -Epnlnf^) , (7) 

n=0 ^-^"^^ 

where p„ represents a "default model" chosen heuristically. The moments 
that enter into the series Figure 1, obtained from simulations if necessary, 
are the information typically used. We then maximize this information en- 
tropy subject to the constraints that the probabilities reproduce the avail- 
able information. The formal maximization of this entropy gives probabil- 
ities 




Pj oc pj exp y- Cfe (^^J j (8) 

where the Cfc are Lagrange multipliers to be adjusted so that the proba- 
bilities finally reproduce the information given initially. The machinery for 
doing this can be developed straightforwardly. For example, the normal- 
ization of the probabilities can be deferred at intermediate stages of the 
calculation. Then the final thermodynamic result can be given in terms of 
the required normalization factor 

This is suggestive of the calculation of a partition function for a modest- 
sized set of states with effective interactions. The interesting questions then 
involve the predictions extracted from the pj for properties other than the 
given information. In our case, the property of first interest is pAjj,. 

We can make the comforting observation that use of only (?7,)o and the 
natural default model pj oc 1/j! produces the Poisson distribution as ex- 
pected; Ci = ~ ln(n)o and /3AiJ, = (n)o as would be found by retaining only 
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the first term shown in Figure 2. Note that evaluation of the second term 
there with the Poisson distribution gives zero as it should [53]; {Sn^)o = (n)o 
for the Poisson distribution. 



10 
9 
8 
7 
6 
5 
4 
3 











1 

^ 














■■, 1 

— t 


1 






t 












— 1 

■, 1 












* « 


> - < 


1 




i, 
i; 




















1 : 
1 'i 

1 ': 




















1 
1 




































' ( 




>---C 


>--< 


5___< 




>---C 




) 



-1.2 



-1.3 



1.4 



1.5 



23456789 
binomial moments 



Figure 3. Convergence of predictions for hydration free energy and information entropy 
with numbers of binomial moments employed; pj oc 1/j!. The moment information was 
obtained from computer simulation of liquid water [58]. A=3.0A. The filled (open) circles 
use the left (right) scale. 



Figure 3 shows how this prediction scheme works-out when the solvent 
is computer simulated liquid water and the solute is a hard sphere of a size 
appropriate for comparison with a Ne atom. The immediate point is that 
a model based upon the two moments that might be obtained from experi- 
ment, (n)o and (n(n — l))o, fortuitously provides the most satisfactory sim- 
ple prediction of po. The Poisson (one moment) model is not satisfactory. 
Inclusion of moments higher than the second is not advantageous unless 
several higher moments arc available. The behavior of the information en- 
tropy f] suggests that the initial two moments do a good job of describing 
the distribution and the subsequent higher moments are 'uninformative.' 

These points are further remarkable because recent analyses [56, 57] 
have underscored the fact that Percus-Yevick approximate integral equa- 
tions can be derived on the assumption that solvent density fluctuations 
are distributed according to a Gaussian probability functional. The results 
of the present investigations, including particularly [56-59], give further in- 
sight and support to those ideas. The PC approximate integral equation 
theory of hydrophobic effects [42], at its inception a Percus-Yevick ana- 
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logue, can be given a similar basis [56]. The PC theory is thereby given a 
better foundation than was available at its genesis. 

Note that present two-moment model is not precisely an assumption of 
a Gaussian probability functional for a density field. The occupancies here 
are required to be nonnegative integers. That this be true of all subvolumes 
of the observation volume is an important restriction. 



TABLE 2. Approximate evaluation of contributions to the series Figure 2 for 
the circumstances of Figure 3, A=3.0A. The final row gives the result obtained 
directly from simulation [58] . The two right-most columns elaborate a direct idea 
for reorganization of this series, motivated by the form of the single-term result 
l3Afi « — ln[l — (71)0]- Though this resummation improves the prediction after 
each term, after five terms the result is still not as successful as the two-moment 
information theory shown in Figure 3. 



j 


(("))o 






_ «.,7a,-i 
- (1-1/j) 


-ln[l-r,] + EiK-r//fe] 


1 


3.77 


3.77 


3.77 






2 


5.75 


1.36 


5.13 


0.72 


5.42 


3 


4.57 


0.75 


5.88 


0.83 


6.29 


4 


2.04 


0.50 


6.38 


0.89 


6.91 


5 


0.51 


0.36 


6.74 


0.90 


7.21 


00 






7.93 







That specifically binomial moments are involved in the series above 
emphasizes the point that the occupancies must be nonnegative integers. 
Further perspective on the convergence issue is obtained by examination 
of successive contributions to the series Figure 2. Numerical results are 
shown in Table 2. The first and second order terms make a significant 
contribution but by themselves are not close to the full answer. The values of 
the additional terms do not establish a rapid convergence to the known full 
answer. However, wc can view the information theory model as a technique 
for reorganization of the series. 

4.1. ENTROPY CONVERGENCE 

When the typical occupation numbers n are large the granularity of the 
distribution p„ is expected to be less significant, at least near the center 
of the distribution and when viewed on a coarse enough scale. In such 
circumstances the predictions of two-moment information theory models 
are not significantly different than those of the PC theory. It is remarkable 
that a simple calculation along these lines gives a convincing explanation of 
the puzzling and contentious issue of "entropy convergence" in hydrophobic 
hydration [60]. 
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The phenomenon to be explained is the following: entropies of transfer 
of non-polar molecules from gas phase or a non-polar solvent into water 
converge at a temperature of about 400 K to approximately zero entropy 
change. Similar behavior was also seen in the microcalorimetry experiments 
on unfolding of several globular proteins. This behavior is insensitive to the 
particular hydrophobic solute molecule. Since the entropy is a temperature 
derivative of a hydration free energy, the convergence temperature identifies 
a region where graphs of hydration free energy versus temperature are ex- 
tremal, in fact, maximal. Below that region the hydrophobic hydration free 
energy increases with temperature but above that region the hydrophobic 
hydration free energy decreases as the temperature is raised. 

The two-moment information theory model above was applied to this 
problem for hard sphere solutes in water with the heuristic modification 
that a flat default model was used; pj oc constant for i<]rnax and zero oth- 
erwise. This latter adjustment was found empirically to give slightly better 
hydration free energies. The results of the model and simulation calcula- 
tions accurately agreed on the temperature dependence of the hydration 
free energies. To analyze this agreement the information theory model was 
simplified to a continuous Gaussian distribution that then gives 



V is the observation volume and psat{T) is the liquid density along the vapor 
saturation curve so that (n)o = Psat{T)v. ((5n^)o was found to be insensi- 
tive to temperature for the relevant conditions. Further, the first term of 
Eq. (11) is larger than the second. Thus the non-monotonic behavior of the 
free energy with temperature and the entropy convergence is a consequence 
of the non-monotonic variation of Tp^atiT)'^ with temperature. The only 
molecular parameter to complicate matters is the volume v and with this 
formula v does not affect the entropy convergence temperature [60]. Thus 
the temperature of entropy convergence is about the same for a wide family 
of solutes. 

The physical point is: the entropy convergence phenomenon occurs for 
water because of the low and temperature insensitive values of {5n?)Q. In 
fact, the isothermal compressibility of water at low pressure has a mini- 
mum value at T=319 K. That temperature differs substantially from the 
observed entropy convergence temperatures but it is not necessary that 
these temperatures be approximately equal, just that ((5n^)o be insensitive 
to temperature in the region d\TpsatiTY\/dT «0. 

The technical point of Eq. (10) is: this formula is simple and effective 
but how it is obtained from the series Figure 2 is not simple. Part of the 
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Tpsat{Tf{kBV^/2{5n\} + T{k^ \n{2TT{5nW/2] 



(11) 
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complication is that two additional twists have been interjected, the flat 
default model and the continuous approximation. 

4.2. PRESSURE DENATURATION OF PROTEINS 

It is a common view, based upon our current understanding, that hydropho- 
bic effects provide a nonspecific, cohesive stabilization of compact protein 
structures. However, it has been argued [62] that our current understand- 
ing of hydrophobic effects is not consistent with the experimental facts of 
pressure denaturation of globular proteins. The information theory model 
of the previous section was applied also to study hydration free energies 
and potentials of mean force (pmfs) for two and three hydrophobic spheri- 
cal solutes in water as a function of pressure [63]. As is well known, those 
pmfs exhibit contact and solvent-separated minima corresponding, respec- 
tively, to cases where the hydrophobic spheres contact each other or where 
a water molecule intervenes. It was found that increasing pressure shifted 
the free energy balance of those two cases towards the solvent-separated 
circumstance. This suggested an intercalation mechanism for pressure de- 
naturation: as the pressure of the liquid is raised, water molecules are forced 
into protein structure. A similar point of view can be taken of the formation 
of clathrate hydrates at elevated pressures: at low pressure hydrophobic ef- 
fects lead to close contacts and to clustering of hydrocarbon gases dissolved 
in water. Pressure increases stabilize the crystalline phase that eliminates 
close solute contacts. If attention is focused on the hydrocarbon material 
this behavior might seem counter-intuitive because the hydrocarbon mate- 
rial seems to expand. But the thermodynamic principle is that increasing 
pressure stabilizes the phase of lower volume. Therefore, we conclude that 
the system may be packed more efficiently and have a lower total volume 
when water molecules are intercalated into the hydrocarbon clusters. These 
topics will surely be the subject of further research. 

5. Concluding Comments 

Identification of some generalizations and future directions for these theo- 
ries will provide concluding comments. Firstly, we note that the generaliza- 
tion of these ideas to treat continuous, rather than only hard core repulsive, 
solute-solvent interactions is known [64]. 

Secondly, we note the importance in the context of the aqueous solu- 
tions of our restriction here to small molecule solutes. It is well recognized 
that treatment of larger solutes requires consideration of the multiphasic 
character of these solutions on large length scales [2, 58, 63]. For large 
enough hard sphere solutes dissolved in water close to phase coexistence, 
the possibility that the solvent will pull away from the solute surface re- 
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quires specific attention. Further subtleties arise when the solute-solvent 
interactions are not just repulsive but include attractive interactions too 
[65, 66]. These issues will surely be the subject of further research in the 
area of hydrophobic effects. It seems likely that an appropriately designed 
default model should be able to describe such effects in a physical manner. 

Thirdly, we note that these approaches provide some unanticipated an- 
swers [64] to questions such as "How is water different from hydrocarbon 
liquids as a solvent for nonpolar solutes?" The importance of the low and 
temperature insensitive values of the isothermal compressibility of liquid 
water is noteworthy. However, such answers are not in the format that is 
most often intended when such questions are asked. Most often such ques- 
tions solicit information about particular patterns of solvent structure in 
the neighborhood of a hydrophobic solutes. Some groundwork has been laid 
for consideration of those detailed structural issues in a format consistent 
with the discussion of this paper [67]. Pursuit of answers about the detailed 
structural issues and their relevance to hydrophobic effects will surely be 
the subject of future research. 

Finally, we note again the relevance of these ideas to the classic problems 
of packing in liquids. The importance of these issues is reflected in the 
significance of the hard sphere fluid system to our understanding of liquids. 
Hard core model systems may not be directly realistic. But it continues to 
surprise that when attention is directed to new physical problems, e.g. the 
thermodynamics and structure of glasses or folded proteins, understanding 
of basic packing problems is again requested. Such problems surface quite 
broadly [66-72]. It would be interesting to see the ordering phase transitions 
associated with packing problems, crystallization and liquid crystallization, 
analyzed on these bases. An initial step along such lines for the hard sphere 
fluid has been taken [75] but more work is deserved. 

We hasten to add that the results so far have not superseded previous 
theoretical results. But this new approach offers the possibility of better, 
more physical understandings of packing problems in the equilibrium sta- 
tistical mechanics of non-crystalline materials and the previous theories of 
them. This approach has achieved new understanding for the problems of 
primitive hydrophobic effects. By exploiting information external to con- 
ventional theories, even simulation data, and by proposing a pattern for 
utilizing that information, these approaches begin to respond to Ander- 
sen's [76] request for a 'theory of theories.' 
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